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q_) , A finite spin system invariant under a symmetry group G is a very illustrative example of the finite 

. group action on a set of mappings f:X—>Y. In the case of spin systems X is a set of spin carriers 

and Y contains 2s + 1 z-components — s < m < s for a given spin number s. Orbits and stabilizers 
are used as additional indices of the symmetry adapted basis. Their mathematical nature does not 
lead to smaller eigenproblems, but they label states in a systematic way. Some combinatorial and 
group-theoretical structures, like double cosets and transitive representations, appear in a natural 
way. In such a system one can construct general formulas for vectors of symmetry adapted basis and 
matrix elements of operators commuting with the action of G in the space of states. Considerations 
presented in this paper should be followed by detailed discussion of different symmetry groups (e.g. 
I 1 cyclic of dihedral ones) and optimal implementation of algorithms. The paradigmatic example, i.e. 

a finite spin system, can be useful in investigation of magnetic macromolecules. 
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I. INTRODUCTION 



Polynuclcar clusters provide magnetic materials with a scale intermediate between that of isolated dimers or trimcrs 
j and that of bulk materials These magnetic materials exhibit features on a mesoscopic scale, so they may show 
quantum effects coexisting with classical behavior. In addition, large assemblies of spins are interesting as real 
objects on which one can test theoretical models with a finite number of spins. Polynuclear magnetic aggregates have 
well denned crystal structures allowing quantitative comparison of theoretical results with experimental ones. 
Unlike other assemblies of small magnetic particles a typical sample of a molecular magnetic compound is composed 
t-H , of nominally identical non-interacting magnets with a unique set of chemically determined parameters. They are 
complex organometallic systems, too difficult to investigate them by ab initio methods, applicable to simple metal 
,— 1 \ cluster. 

The aim of this paper is to present and discuss combinatorial structure of finite spin models. Such investigation lead 
to general analytical formulas for the irreducible basis and matrix elements of any operator commuting with all sym- 
metry operators P{g), g £ G, The most important is a Hamiltonian of the system in question, but also other operators 
' can be taken into account. Since S 2 commutes with the isotropic Heisenberg Hamiltonian, then the total spin S can be 

used as an additional label of states. Finite spin models have been used in condensed matter physics for many years 
■ and have been applied to various models of magnetic materials: from the one-dimensional ferromagnetic Ising model 
to the three-dimensional Heisenberg antiferromagnets. Recent developments and discoveries have stimulated interests 
in finite spin models, since they can be applied to meso- and nanoscopic systems, especially magnetic macromolecules 
with well-determined symmetry groups like [NaFeefOCHgji^CiTHisO^elClCU (Fe 6 ), [FeioCOCHs^ot^^OaCTjxo] 
(Feio), Mni 2 acetate or [Mn(hfac) 2 NITPh] 6 (Mn 6 ) || g, § §. Considering macromolecules with many magnetic cen- 
ters and large spin numbers (up to s — 5/2) one deals with a very large number of states even for a small number 
£j , of magnetic ions N. Moreover, this number increases rapidly with increasing N and s. However, group-theoretical 
and combinatorial methods yield decomposition of a space of states into subspaces with relatively small dimensions. 
Therefore, magnetic properties of macromolecules can be investigated very deeply and model parameters can be fitted 
in more accurate way. This paper, based on the results of some earlier papers jR h], is one of the steps necessary in 
the study considerations of small (nanoscopic) rings with relatively large spin number. 

It should be also stressed that the approach described and solutions suggested can be used in any problem where 
a Hamiltonian symmetry group G (or, at least, its subgroup) permutes basis vectors. It means that a permutation 
representation P of G is determined or, in other words, a group action of G on a basis B of quantum space of states 
is introduced. This immediately gives rise to such structures and concepts as transitive representations, orbits, stabi- 
lizers, double cosets etc. Permutation representations themselves constitute a special case of induced representation 
and are applied in many branches of physics |9| . Some problems related to permutation and transitive representations 
were studied by Lulek et al. P, [Hi and their results, in various extent, are used here. 
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II. FINITE SPIN MODELS AND SIMPLE COMBINATORIAL OBJECTS 

A system of iV spin s gives rise to (2s + 1) Ising configurations 

\mi m 2 . . . tun) , rrij S {— s, — s + 1, . . . , s} . 
They can be considered as mappings 

fi:{l,2, ...,N}-> {-s,-s+ l,...,s}:j ^ fi(j) = m . 

Magnetization of the system, M = X)j=i MO) — S^=i m jj ^ s an example of an additive weight. Therefore, all 
combinatorial methods of enumeration with weights can be applied [p"l| . Moreover, each mapping (configuration) is 
also characterized by a content, i.e. by numbers 

[k] := [fc_ s , . . . , k s ] 

being numbers of solutions of equations 

= m; me {— s, — s + 1, . . . , s} . 

In fact [k] is a non-ordered partition of N. A number of configurations with the same content is given by a polynomial 
coefficient 

AH (V s k m Y- 

n(\k\) — ' — KL ^ m =- s m > 

Hence, a number of configurations (mappings) with a given weight (magnetization) is 

n(N,M)= Yl "CM)" C 1 ) 

Let 

£) s (g) £> s (g) . . . £^ = n(S)D s . 

N times 5=0 or V2 

be a decomposition of the iVth tensor power of the (spin) irrep D s . Then 

n(S = Ns) = n(N, M — Ns) = 1 ; n(S) = n(N, M = S) — n{N, M = S + 1) for S < Ns . 

Some recurrences for n([fc]) and n(M) can derived from those for polynomial coefficients and depicted by analogs 
of the Pascal triangle. In Fig. [l]a the recurrence 

s 

n ([k}) = X! n (i k -s,k-s+i, ■ ■ ■ ,kj - 1, . . .,&„]) 

is illustrated by arrows for n[l, 1, 1] = n[0, 1, 1] + n[l, 0, 1] + n[l, 1, 0]. Introducing this relation to Eq. (Q) one obtains 
§: 

S 

n(N,M)= Y n(N -l,M-m). 

771 — — S 

It is shown by arrows in Fig. lb in the case N = 3, M = -1: n(3, -1) = u(2, -2) + n(2, -1) + n(2, 0). The dotted 
lines in Fig. la indicate polynomial coefficients summed in Eq. (|l|). After summation one obtains numbers presented 
in the last row in Fig. lb. Formal equations for the number of configurations with given iV and M can be derived by 
methods presented in Kerber's monograph [ fTl) . 
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FIG. 1: Recurrences for a) polynomial coefficients and b) numbers n(M) for s = 1 



III. FINITE GROUP ACTION 

For finite both a group G and a set A an action qX of G on X is defined as a mapping (an external multiplication 
rule in X) 

(G x X) — y X: (g,x) i — > 52; with {9i92)x = 51(52)2; and lga; = a; . 
A stabilizer G x of x G A is a subgroup of G such that 

Gx '■= {5 6 G I 52; = . 
An orbit G(x) of x S A is a subset of A such that 

G(x) := { 3 x I 5 G G} ; |G(x)||G^| = |G| . 
The action gl can be raised to an action on the set Y x of all mappings /: A — > Y (|Y X | = |Y|I X I for finite sets): 
(Gx F X )->Y*: /->/' = /os- 1 , i.e. /'(.t) = /(g-^) . 
Enumeration of mappings needs another subset of A — the so-called fixed points: 

X g = {x G A I gx = x} . 

This set is used, for example, to determine number of mappings with given weight and stabilizer by means of the 
so-called Polya substitution 11 . 



IV. TRANSITIVE REPRESENTATIONS, DOUBLE COSETS AND SO ON 

Any action of G on an orbit G{x) is similar[^U to the action of G on the set of left cosets 

G/G x = {g r G x I \<r < \G(x)\} , 
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where r labels representatives g r of left cosets; this action is defined as 

g(g r G x ) = (gg r )G x . 

If subgroups G x and G y are conjugated in G, then actions g{G/G x ) and g(G/G v ) are similar, ft means that 
nonequivalent actions are labeled by subgroups U representing classes of conjugated subgroups. These actions, 
g(G/U), can be considered as building blocks of any action gX. 

Let £> = {|r) | 1 < r < |G/J7|} be a basis of a (formal) unitary space L constructed as a linear closure (over the 
field of complex numbers) of B. The vectors in B are in a one-to-one correspondence with elements of an orbit G(U), 
i.e. with left cosets g r U . Moreover, a representation P of G in the space L can be defined in a natural way 12 : 



P(g)\r) = \s) if gg r eg s U. 

This representation is denoted as R G:U and called transitive. 

A transitive representation R G ' U is, in a general case, reducible and can be decomposed into a direct sum of 
irreducible representations (irreps) of G: 

rG:U = 0n(I\[/)r ; 0<n(T,U)< [T], 

where [r] denotes dimension of T. Any action qX determines a permutation representation of P in a space spanned 
over vectors \xr), where x labels orbits and r labels representatives of cosets g r G x . Since each finite action can be 
decomposed into action on orbits, then a permutation representation is decomposed into transitive representations. 
A double coset is determined by two subgroups of G and is defined as 

VgU = {vgu \ v £ V : u E U} . 

Let us consider an action of G on the Cartesian product G(y) x G(x) of two orbits with the stabilizers G x = U and 
G(y) = V. This action is defined as 

(9, (y,x)) >-> (gy,gx) 

and has the following properties (x and y are fixed): 

• Each orbit contains an element of the form (y,gx); 

• The stabilizer of (y,gx) is G y n gG x g~ 1- , 



The mapping G(y, gx) 1— > G y gG x is a bijection, so orbits of the introduced action can be labeled by double coset 
representatives. 



V. EQUATIONS FOR MATRIX ELEMENTS 

It can be shown [ p"3[ |l4| that any operator Ji commuting with all operators P{g), J £ G, can be written in a 
quasi-diagonal form with blocks labeled by irreps T. In each block rows (columns) are indexed by a stability group 
U (V), representatives x (y) of orbits with a stabilizer U (V, respectively), and by indices of vectors |r) (|s)) of the 
irrep T. A general form of matrix elements is following: 

\V\G/U\ 

/^,v^(r) = (|(7||F|r 1/2 \Vg d U\{Uxg d \H\Vye G )B^ r {g d ), (2) 

d=l 

where B^ T (gd) depends on the so-called reduction coefficients and can be treated as group-theoretical parameters 
of a model under considerations. Analytic formulas for these coefficients can be determined in the case of cyclic or 
dihedral group G |^|, This equation simplifies for the unit representation To'. 

nm ,l/2\G/U\ 

h Ux , Vy (r ) = ^jyjj J2 (Ux gj \H\Vye G ) , 
If there is only one double coset, represented by ec, then all left coset representatives belong to it and 

W,(To)= (|[^[) ' ( Uxe G\n\Vye G ), 



so it is enough to calculate matrix elements for orbit representatives. 
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TABLE I: Multiplicities n(T) for a ring of N = 12 spins s = 1 with G = D12 and M = 

T Ai A 2 Bi B 2 £1 £2 £3 £4 £5 
n(r) 3179 2987 3107 3056 6136 6158 6140 6160 6136 

VI. SOME REMARKS ON NUMERICAL PROBLEMS 

Application of the formula presented above is not immediate and easy. Of course, in very simple cases, e.g. N = 4 and 
s = 1, one can calculate all matrix elements by hand. Larger systems require proper construction and implementation 
of algorithms. Moreover, one can use additional quantum numbers (good quantum numbers as the magnetization 
and the total spin in the case of isotropic Heisenberg Hamiltonians) to obtain more detailed classification of states. 
Secondly, the labeling scheme of orbit representatives may include partitions of N, which, in fact, label orbits of the 
symmetric group D G. The action of this group does not commute with the Hamiltonian, in a general case, but 
commutes with the operator S z — J>2j=i s j °f the total magnetization. Therefore, if one wants to use the quantum 
number M as an additional index, then it is natural to include partitions in the labeling scheme. 

One of advantages of the formula presented is a fact, that in some cases (especially for one-dimensional irreps T) 
obtained matrix elements are integers. Therefore, considering the eigenproblems for S 2 , which has integer eigenvalues 
(for N odd and half-integer s one can consider 4S 2 instead) , we obtain a set of homogeneous linear equations with 
integer coefficients. This set can be solved exactly, but we have to cope with very large numbers. This can be done 
by means of any multiple precision library, for example GMP ]l7| . Prototypes provided by this library can be also 
used to store configurations in compact (packed) form. Moreover, there are as well procedures for high-precision float 
calculations, so the Hamiltonian eigenproblem can be solved, too. More detailed discussion of numerical problems 
can be found in other papers Jl3|, |i~§|| . 

VII. EXAMPLES 

The proposed methods are especially useful in the case of ring-shaped magnetic molecules with cyclic of dihedral 
symmetry group. Below some results for two molecules are presented: (i) Nii2(O2CMe)i2(chp)i2(H 2 0)6(THF) 6 re- 
ferred to as Ni 12 and (ii) Fe6- The first one has the ferromagnetic ground state wit the total spin S = 12 0, whereas 
the second is an example of a molecular antiferromagnet with S = in the ground state || . 

The M12 cluster contains a ring of twelve s = 1 spin carriers. There are 3 12 = 531441 Ising configurations. Assuming 
isotropic interactions each subspace with a given magnetization — 12 < M < 12 can be considered separately. The 
largest dimension, 73789, has the subspace with M — 0. Ising configurations spanning this space are labeled by the 
following non-ordered partitions of 12 into no more than three parts (number of configurations for each partition is 
given in the parenthesis): [0,12,0] (1), [1,10,1] (132), [2,8,2] (2970), [3,6,3] (18480), [4,4,4] (34650), [5,2,5] (16632), and 
[6,0,6] (924). In the next step orbit of the symmetric group £12 are decomposed into orbits of the dihedral group Di 2 , 
what leads to transitive representations; they are decomposed into irreps of D12. For each irrep one can construct a 
block of the Hamiltonian matrix according with Eq. (||) . The dimensions of these blocks are collected in Table Q; the 
largest one is 6160 for Y = E4.. All eigenvalues (for too) of the Hamiltonian 

12 

H = - y^(JsjS j+ i + guBBs*) 
j'=i 

can be calculated with high numerical precision and thermodynamic properties can be determined. The best fit to 
experimental results presented by Blake et al. [[l9| is achieved for the exchange integral J/ks = 8.5 ± 0.5 K and the 
gyromagnetic ratio g = 2.23 ± 0.01 Q. 

The antiferromagnetic ring Fe6 contains six magnetic ions with s = 5/2, so there are 46656 Ising configurations. 
Numbers n(S) and n(6, M) for M > are presented in Table [□]. Considering the ground state one has to take into 
account 32 non-ordered partitions. Decomposition into orbits of the dihedral group involves only three stabilizers: Ci, 
Di, and D3. Taking into account decompositions of transitive representations into irreps one obtains results collected 
in Table |l]. 

Our aim is to calculate spin correlation in the ground state assuming isotropic Heisenberg interactions. In the 
considered case (N — 6, s — 5/2) the ground state has symmetry V — B\, so the dimension of the eigenproblem 
is 385. However, we solve it in two steps: (i) at first we determine and orthonormalize eigenvectors of S 2 operator 
with the eigenvalue zero; (ii) next we transform the Hamiltonian matrix to obtain a block labeled by Y — B\ and 
M = S = 0. The eigenproblem for such small matrix (20 x 20) can be easily solved with high numerical precision. We 
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TABLE II: Numbers n(S) and n(6, M) for S, M = 0, 1, . . . , 15 
g or M 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

n(S) 111 315 475 575 609 581 505 405 300 204 126 70 35 15 5 1 

(2S + l)n(S) 111 945 2375 4025 5481 6391 6565 6075 5100 3876 2646 1610 875 405 145 31 
n(6,M) 4332 4221 3906 3431 2856 2247 1666 1161 756 456 252 126 56 21 6 1 

TABLE III: Numbers n(Y) for M = 0, 1 and S = 0; the last one is simply calculated as the difference of the previous two 

T Ai A 2 Bi B 2 £i E 2 

n(T), M = 385 339 385 339 721 721 
n(r), M = 1 383 325 365 334 699 708 
n(r), S = 2 14 20 5 22 13 

determine not only eigenvalues, but eigenstates, too, what allows to calculate spin correlations. The first three low- 
lying states have energies per spin (in the unit of J): —7.322, —5.761, —4.831, what agrees with the results presented 
by Lascialfari et al. Q. The spin correlations per bond uj z (r) — {J2j sjs| +r }/6 in these states are following: 

r = 1 r = 2 r=3 

-2.440817 +1.958689 -1.952410 

-1.920315 +1.056202 -1.188439 

-1.610233 +0.226289 -0.148779 

Note that absolute value of the nearest-neighbor z-z correlation in the ground state is smaller than in the Neel state 
\N) = \s, — s, s, — s, s, — s), when one obtains the value of 6.25. However, the total spin correlation, equal 3uj z (1), 
is stronger than in the Neel state, since x-x and y-y correlations are equal to zero in this state. This fact shows 
importance of other sates with AI — and follows from quantum character of the model considered. 

VIII. FINAL REMARKS 

In this report we present briefly a way from (finite) group action to calculations of spin correlations. The group- 
theoretical and combinatorial object appearing in this course help us to label Ising configurations in the systematic 
way and to obtain relatively simple form of the matrix elements. To exploit all possible simplifications we have 
assume isotropic interactions. However, there is no limit imposed on its range: one can easily include next-nearest 
neighbors or biquadratic terms (sj-Sfc) 2 . In the case of uniaxial anisotropy the total spin S is not a good quantum 
number, but even in this case thermodynamic properties of small magnetic macromolecules can be determined (the 
total spin S was not used to obtain results presented in Q). Such model is frequently assumed to achieve a better 
fit to experimental results, for example in the case of Mni2 acetate molecule Q. The XYZ models, when the total 
magnetization is not a good number, are most difficult. However, the Ising configurations still can be labeled by 
partitions and stabilizers, but the Hamiltonian matrix cannot be decomposed into blocks labeled by M; of course 
the label T can be used as long as the Hamiltonian commutes with all operators P{g). It should be also stressed 
that some results presented are obtained using procedures included in the GMP package. For example, we are able to 
solve exactly a set of homogeneous linear equation with integer coefficients. The methods proposed has been applied 
to magnetic macromolecules, however calculations in the thermodynamic limit are also very important, since there 
are still some open problems related with linear magnets {e.g. critical exponents near T = OK in the case of ferro- 
magnetic ordering |Bl[ ). 
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